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Abstract 


Traditionally,  radar  signals  have  been  analysed  in  either  the  time  or  the  frequency 
domain.  Joint  time-frequency  representations  characterize  signals  over  a 
time-frequency  plane.  They  thus  combine  time-domain  and  frequency-domain  analyses 
to  yield  a  potentially  more  revealing  picture  of  the  temporal  localization  of  a  signal’s 
spectral  components.  Therefore,  for  air  target  returns  with  time-varying  frequency 
content,  the  joint  time-frequency  representations  offer  a  powerful  analysis  tool.  A 
concise  review  of  time-frequency  transforms  is  provided  as  background  and  is  needed 
to  appreciate  how  time-frequency  processing  methods  can  improve  conventional  time 
or  frequency  processing  methods.  The  report  then  describes  and  illustrates  the 
advantages  of  using  joint  time-frequency  techniques  to  analyze  a  multi-component 
signal,  a  noisy  signal,  and  experimental  aircraft  data.  Finally,  we  use  time-frequency 
analysis  techniques  for  the  detection  of  maneuvering  aircraft  using  HF  radar  in  heavily 
cluttered  regions.  We  compare  the  ability  of  different  time-frequency  transforms  to 
resolve  several  experimental  aircraft  returns.  The  results  clearly  demonstrate  that 
time-frequency  analysis  techniques  can  significantly  improve  the  detection 
performance  of  the  HF  radar  and  add  considerable  physical  insight  over  what  can  be 
achieved  by  conventional  Fourier  transform  methods  currently  used  by  HF  radars. 


DRDC  Ottawa  TM  2003-090 


I 


Resume 


Par  le  passd,  les  signaux  radar  ont  ete  analyses  soit  dans  le  domaine  temporel,  soit  dans 
le  domaine  frequentiel.  Les  representations  temps-frequence  mixtes  montrent  les 
caracteristiques  des  signaux  sur  un  plan  temps-frequence.  Elies  combinent  ainsi  les 
analyses  du  domaine  temporel  et  du  domaine  frequentiel  pour  donner  une  image 
potentiellement  plus  revelatrice  de  l’emplacement  temporel  des  composantes  spectrales 
d’un  signal.  Par  consequent,  les  representations  temps-frequence  mixtes  constituent  un 
outil  d’analyse  puissant  des  signaux  dont  le  contenu  frequentiel  varie  dans  le  temps. 
Nous  presentons  une  breve  etude  des  transformations  temps-frequence  afin  de  foumir  la 
base  necessaire  pour  apprecier  la  fa9on  dont  les  methodes  de  traitement 
temps-frequence  permettent  d’ameliorer  les  methodes  classiques  de  traitement  temporel 
ou  frequentiel.  Ensuite,  le  rapport  decrit  et  illustre  les  avantages  de  l’utilisation  des 
techniques  temps-frequence  mixtes  pour  1’ analyse  d’un  signal  a  plusieurs  composantes, 
d’un  signal  perturbe  par  le  bruit  et  de  donnees  experimentales  d’aeronef.  Enfin,  nous 
utilisons  les  techniques  d’analyse  temps-frequence  pour  la  detection  d’un  aeronef 
effectuant  des  manoeuvres  en  ayant  recours  au  radar  HF  dans  des  zones  a  niveau  de 
clutter  el  eve.  Nous  comparons  la  capacite  de  differentes  transformations 
temps-frequence  pour  resoudre  plusieurs  signaux  experimentaux  d’aeronef.  Les 
resultats  montrent  clairement  que  les  techniques  d’analyse  temps-frequence  permettent 
d’amdliorer  de  fa?on  appreciable  l’efficacite  de  detection  du  radar  HF  et  d’acquerir 
beaucoup  d’autres  connaissances  physiques  par  comparaison  aux  methodes  a 
transformation  de  Fourier  classiques  qui  sont  actuellement  utilisees  par  les  radars  HF. 
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Executive  summary 


One  of  the  central  problems  in  exploiting  High  Frequency  (HF)  radar  data  is  the 
analysis  of  a  time  series.  The  problem  at  hand  is  how  to  extract  the  information  present 
in  the  data  and  use  it  to  its  full  potential.  In  order  to  address  this  problem  we  turn  to  the 
field  of  signal  analysis  and  data  representations.  Traditionally,  radar  signals  have  been 
analysed  in  either  the  time  or  the  frequency  domain.  The  Fourier  Transform  is  at  the 
heart  of  a  wide  range  of  techniques  that  are  used  in  HF  radar  data  analysis  and 
processing.  Mapping  the  data  into  the  temporal  frequency  domain  is  an  effective  way  of 
recording  the  data  such  that  their  global  characteristics  can  be  assessed.  However,  the 
change  of  frequency  content  with  time  is  one  of  the  main  features  we  observe  in  HF 
radar  data.  Because  of  this  change  of  frequency  content  with  time,  radar  signals  belong 
to  the  class  of  non-stationary  signals.  The  analysis  of  non-stationary  signals  requires 
technique  that  extend  the  notion  of  a  global  frequency  spectrum  to  a  local  frequency 
description.  The  spectral  energy  density  function  that  is  obtained  by  means  of  a  Fourier 
Transform,  the  so-called  power  spectrum,  shows  the  frequencies  that  are  present  in  our 
data,  but  does  not  reveal  where  changes  in  the  frequency  content  occur.  Consequently, 
for  the  interpretation  of  radar  data  in  terms  of  a  changing  frequency  content,  we  need  a 
representation  of  our  data  as  a  function  of  both  time  and  frequency.  Only,  quite 
recently,  the  joint  time-frequency  representation  of  signals  has  become  a  major  area  of 
research  in  radar  signal  processing. 

The  time-frequency  representations  characterize  signals  over  a  time-frequency  plane. 
They  thus  combine  time-domain  and  frequency-domain  analyses  to  yield  a  potentially 
more  revealing  picture  of  the  temporal  localization  of  a  signal’s  spectral  components. 
They  may  also  serve  as  a  basis  for  signal  detection,  imaging,  characterization,  coding, 
and  processing.  A  complete  and  comprehensive  theory  for  joint  time-frequency 
analysis  does  not  yet  exist.  There  is  no  unique  time-frequency  representation  of  a  signal 
that  satisfies  all  the  properties  of  a  physically  correct  joint  time-frequency  energy 
density  function.  However,  discarding  the  requirement  that  all  properties  must  be 
satisfied  in  one  time-frequency  representation,  a  class  of  joint  time-frequency 
representations  can  be  derived  that  serves  as  a  model  of  a  local  power  spectrum. 

In  this  report,  a  concise  review  of  time-frequency  transforms  is  provided  as  background 
and  is  needed  to  appreciate  how  time-frequency  processing  methods  can  improve 
conventional  time  or  frequency  processing  methods.  The  report  then  describes  and 
illustrates  the  advantages  of  using  joint  time-frequency  techniques  to  analyze  a 
multi-component  signal,  a  noisy  signal,  and  experimental  aircraft  data.  Finally,  we  use 
time-frequency  analysis  techniques  for  the  detection  of  maneuvering  aircraft  using  HF 
radar  in  heavily  cluttered  regions.  We  compare  the  ability  of  different  time-frequency 
transforms  to  resolve  several  experimental  aircraft  signals.  The  results  clearly 
demonstrate  that  time-frequency  analysis  techniques  can  significantly  improve  the 
detection  performance  of  the  HF  radar  and  add  considerable  physical  insight  over  what 
can  be  achieved  by  conventional  Fourier  transform  methods  currently  used  by  HF 
radars. 
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Sommaire 


L’ analyse  d’une  serie  chronologique  constitue  l’un  des  principaux  problemes  eprouves 
avec  les  donnees  du  radar  haute  frequence  (HF).  Le  probleme  qui  se  pose  consiste  a 
trouver  un  moyen  pour  extraire  1 ’information  contenue  dans  les  donnees  et  pour 
l’utiliser  a  son  plein  potentiel.  Pour  resoudre  ce  probleme,  nous  faisons  appel  au 
domaine  de  l’analyse  des  signaux  et  de  la  representation  des  donnees.  Par  le  passe,  les 
signaux  radar  ont  ete  analyses  soit  dans  le  domaine  temporel,  soit  dans  le  domaine 
frequentiel.  La  transformation  de  Fourier  est  au  coeur  d’une  vaste  gamme  de  techniques 
qui  sont  utilisees  pour  Panalyse  et  le  traitement  des  donnees  du  radar  HF.  La  mise  en 
correspondance  des  donnees  dans  le  domaine  des  frequences  temporelles  constitue  un 
moyen  efficace  pour  enregistrer  les  donnees  de  telle  fafon  que  leurs  caracteristiques 
globales  puissent  etre  evaluees.  Cependant,  la  variation  du  contenu  frequentiel  avec  le 
temps  est  une  des  principales  caracteristiques  que  nous  observons  dans  les  donnees  du 
radar  HF.  En  raison  de  cette  variation  du  contenu  frequentiel  avec  le  temps,  les  signaux 
radar  appartiennent  a  la  categorie  des  signaux  non  stationnaires.  Pour  l’analyse  des 
signaux  non  stationnaires,  il  faut  avoir  recours  a  des  techniques  qui  elargissent  la  notion 
d’un  spectre  de  frequences  global  a  une  description  de  frequence  locale.  La  fonction  de 
densite  d’energie  spectrale  obtenue  au  moyen  d’une  transformation  de  Fourier,  et 
appelee  le  spectre  de  puissance,  montre  quelles  frequences  sont  contenues  dans  nos 
donnees,  mais  n’indique  pas  ou  se  produisent  les  changements  du  contenu  frequentiel. 
Par  consequent,  pour  1’ interpretation  des  donnees  radar  du  point  de  vue  d’un  contenu 
frequentiel  changeant,  nous  avons  besoin  d’une  representation  de  nos  donnees  en 
fonction  du  temps  et  de  la  frequence.  Ce  n’est  que  tout  recemment  que  la  representation 
temps-frequence  mixte  des  signaux  est  devenue  un  sujet  important  de  recherche  dans  le 
traitement  des  signaux  radar. 

Les  representations  temps-frequence  mixtes  montrent  les  caracteristiques  des  signaux 
sur  un  plan  temps-frequence.  Elies  combinent  ainsi  les  analyses  du  domaine  temporel 
et  du  domaine  frequentiel  pour  donner  une  image  potentiellement  plus  revelatrice  de 
l’emplacement  temporel  des  composantes  spectrales  d’un  signal.  Elies  peuvent  aussi 
servir  de  base  pour  la  detection  des  signaux,  l’imagerie,  la  caracterisation,  le  codage  et 
le  traitement.  II  n’existe  pas  encore  de  theorie  complete  et  detaillee  pour  l’analyse 
temps-frequence  mixte.  II  n’existe  pas  de  representation  temps-frequence  unique  d’un 
signal  qui  satisfait  a  toutes  les  proprietes  d’une  fonction  de  densite  d’energie 
temps-frequence  mixte  correcte  du  point  de  vue  physique.  Cependant,  si  on  fait 
abstraction  de  l’exigence  selon  laquelle  il  faut  satisfaire  a  toutes  les  proprietes  dans  une 
seule  representation  temps-frequence,  on  peut  etablir  une  classe  de  representations 
temps-frequence  mixtes  qui  sert  de  modele  pour  un  spectre  de  puissance  local. 

Dans  ce  rapport,  nous  presentons  une  breve  etude  des  transformations  temps-frequence 
afin  de  foumir  la  base  necessaire  pour  apprecier  la  fafon  dont  les  methodes  de 
traitement  temps-frequence  permettent  d’ameliorer  les  methodes  classiques  de 
traitement  temporel  ou  frequentiel.  Ensuite,  le  rapport  decrit  et  illustre  les  avantages  de 
l’utilisation  des  techniques  temps-frequence  mixtes  pour  l’analyse  d’un  signal  a 
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plusieurs  composantes,  d’un  signal  perturbe  par  le  bruit  et  de  donnees  experimentales 
d’aeronef.  Enfin,  nous  utilisons  les  techniques  d’analyse  temps-frequence  pour  la 
detection  d’un  aeronef  effectuant  des  manoeuvres  en  ayant  recours  au  radar  HF  dans  des 
zones  a  niveau  de  clutter  eleve.  Nous  comparons  la  capacite  de  differentes 
transformations  temps-frequence  pour  resoudre  plusieurs  signaux  experimentaux 
d’aeronef.  Les  resultats  montrent  clairement  que  les  techniques  d’analyse 
temps-frequence  permettent  d’ameliorer  de  fa?on  appreciable  l’efficacite  de  detection 
du  radar  HF  et  d’acquerir  beaucoup  d’autres  connaissances  physiques  par  comparaison 
aux  methodes  a  transformation  de  Fourier  classiques  qui  sont  actuellement  utilisees  par 
les  radars  HF. 
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1.  Introduction 


One  fundamental  mathematical  tool  employed  in  radar  signal  processing  is  a  transform. 
When  asked  to  multiply  the  Roman  numerals  LXIV  and  XXXII,  only  a  few  of  us  will 
be  able  to  give  the  correct  answer  right  away.  However,  if  the  Roman  numerals  are  first 
translated  into  Arabic  numerals,  64  and  32,  then  all  of  us  can  get  2048  immediately. 

The  process  of  converting  the  unfamiliar  Roman  numerals  into  common  Arabic 
numerals  is  a  typical  example  of  transforms  [1].  By  properly  applying  transforms,  we 
can  simplify  calculations  or  make  certain  attributes  of  the  signal  explicit. 

One  of  the  most  popular  transforms  known  to  scientists  and  engineers  is  the  Fourier 
transform  that  converts  a  signal  from  the  time-domain  to  the  frequency-domain.  In  fact, 
the  Fourier  transform  is  not  simply  a  mathematical  tool  to  make  calculations  easier.  It 
also  acts  as  a  mathematical  prism  to  break  down  a  signal  into  a  group  of  waveforms 
(different  frequencies),  as  a  prism  breaks  up  light  into  a  color  spectrum.  The  Fourier 
transform  is  so  powerful  that  people  tend  to  apply  it  everywhere  without  noticing  one 
fundamental  difference  between  the  mathematical  prism  and  a  real  prism.  The 
spectrum  produced  by  the  prism  in  the  morning  is  different  from  that  in  the  evening. 

We  may  say  that  the  prism  gives  instantaneous  spectra.  Using  a  prism  to  examine 
spectra  of  light,  there  is  no  need  for  the  information  about  light  that  existed  million 
years  ago  and  the  light  that  will  be  there  tomorrow.  However,  this  is  not  the  case  for 
Fourier  transform.  To  compute  the  Fourier  transform,  we  not  only  need  previous 
information,  but  also  information  that  have  not  yet  occurred.  The  spectrum  computed 
by  the  Fourier  transform  is  the  spectrum  averaged  over  an  infinitely  long  time  before 
the  present  to  an  infinitely  long  time  after  the  present!  [2] 

The  time  series  signals  of  radar  targets  have  been  traditionally  analyzed  in  the 
frequency-domain.  The  transformation  of  a  time  series  into  a  frequency-domain  using 
the  Fourier  transform  provides  useful  information  regarding  the  target.  However, 
Fourier  spectral  analysis  is  not  without  its  shortcomings.  The  most  prominent 
shortcoming  of  Fourier  spectral  analysis  is  its  inability  to  properly  represent  frequency 
contents  that  changes  over  an  observation  period.  'The  most  widely  used  signal 
processing  tool  is  the  FFT  (fast  Fourier  transform);  the  most  widely  misused  signal 
processing  tool  is  also  the  FFT"  [2],  Fourier  transform-based  techniques  are  effective 
as  long  as  the  frequency  contents  of  the  signal  do  not  change  with  time.  However,  when 
the  frequency  contents  of  the  data  samples  evolve  over  an  observation  period,  other 
alternative  methods  should  be  considered.  Joint  time-frequency  representations  were 
developed  to  overcome  this  limitation  of  Fourier  spectral  analysis  [3-4], 

The  joint  time-frequency  analysis  has  been  a  topic  of  much  interest  in  the  signal 
processing  community  in  the  past  decade  .  Over  the  past  ten  years,  radar  researchers 
have  also  investigated  time-frequency  transforms  as  a  unique  tool  for  radar-specific 
signal  analysis  and  image  processing  applications  [2,5-9],  Both  traditional 
time-frequency  techniques,  as  well  as  the  new  tools  developed  in  the  signal  processing 
community,  have  been  applied  to  various  radar  problems.  Like  the  developments  in 
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other  fields,  such  as  underwater  acoustics  and  speech  processing,  it  was  found  that 
time-frequency  transforms  provide  additional  insight  into  the  analysis,  interpretation, 
and  processing  of  radar  signals  that  is  sometimes  superior  to  what  is  achievable  in  the 
traditional  time  or  frequency  domain  alone.  The  specific  applications  where 
time-frequency  transforms  have  been  used  include  signal  analysis  and  feature 
extraction,  motion  compensation  and  image  formation,  signal  denoising,  imaging  of 
moving  targets,  and  detection  of  moving  targets  [2,5-9], 

Joint  time-frequency  representations  transform  a  one-dimensional  time  signal  into 
two-dimensional  representation  in  the  time-frequency  plane.  In  this  representation, 
time  varying  frequency  can  be  easily  displayed  and  studied.  An  important  advantage  of 
time-frequency  representations  is  the  ease  of  identifying  target  signals.  A  waveform 
with  multiple  signals  or  noise  can  become  too  complicated  for  visual  analysis  of  its 
waveform.  Even  the  Fourier  spectral  analysis  can  become  garbled  beyond  recognition. 
Since  time-frequency  representations  display  time- varying  frequencies,  the  signals  are 
distinguishable  from  die  noise.  The  nature  of  each  signal  can  be  easily  seen. 

This  report  demonstrates  the  usefulness  and  effectiveness  of  time-frequency 
representations  and  compares  various  time-frequency  transforms.  Three  cases  are 
examined  to  show  the  usefulness  of  time-frequency  representations.  The  first  case 
involves  chirp  signals  that  cannot  be  identified  using  Fourier  spectral  analysis.  The 
second  case  involves  a  signal  with  noise.  The  data  in  these  cases  are  simulated  with 
complex  values  to  be  similar  to  real  data.  The  third  case  shows  the  overall  effectiveness 
of  time-frequency  representations  in  signal  analysis.  This  case  uses  real  aircraft  data. 

There  are  many  different  time-frequency  representations  available  in  the  literature  for 
analysis  [2,5-10],  Each  transform  offers  advantages  and  disadvantages.  Certain 
transforms  provide  excellent  resolution,  but  are  slow  to  calculate.  Others  provide  poor 
resolution,  but  are  quick  to  calculate.  Some  transforms  provide  a  good  balance  of  speed 
and  resolution,  but  have  interference  cross-terms.  Although  there  are  many  different 
transforms  available,  we  compare  the  ability  of  only  twelve  transforms  to  resolve 
several  real  aircraft  signals.  The  main  objective  of  this  report  is  to  demonstrate  the 
effectiveness  of  time-frequency  transforms  for  detecting  maneuvering  air  targets  in 
heavily  cluttered  environment.  Twelve  transforms  that  are  considered  for  comparison 
and  evaluation  are  the  adaptive  energy  distribution,  Bom- Jordan  distribution,  binomial 
distribution,  Choi-Williams  distribution,  Gabor  representation,  Margenau-Hill 
spectrogram,  smoothed  Wigner-Ville  distribution,  spectrogram  distribution, 
Zhao-Atlas-Marks  distribution,  reassigned  Gabor  representation,  reassigned  Morlet 
scalogram,  and  the  reassigned  spectrogram  distribution.  Other  time-frequency 
transforms  that  are  available  in  the  literature  will  be  evaluated  and  compared  in  a  future 
report.  These  twelve  transforms  are  compared  with  experimental  radar  data.  There  are 
several  examples  of  target  returns  with  varying  behavior  to  test  the  transforms  in 
different  situations.  Each  transform  is  then  evaluated  individually  and  against  the  group. 
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2.  Advantages  of  Time-Frequency  Representations 


Time-frequency  representations  can  drastically  simplify  the  interpretation  of  signals. 
These  representations  provide  a  much  clearer  image  of  the  signals  due  to  their  ability  to 
display  frequency  that  changes  over  time.  To  show  this  advantage  more  clearly,  we  will 
examine  three  cases.  The  first  case  consists  of  multiple  chirp  signals.  The  second  case 
consists  of  a  signal  present  in  noise.  These  cases  use  simulated  data  with  complex 
values.  The  third  case  uses  experimental  aircraft  data  to  show  the  usefulness  of 
time-frequency  representations  for  signal  analysis. 

2.1  Detecting  Multiple  Signals 


Real  Part  of  Signal  in  Time 


Figure  1:  Time  series  of  multi-component  signal 

There  are  multiple  chirps  present  in  the  signal  whose  time  series  is  displayed  in  Figure 
1 .  The  spectrum  of  the  signal  can  be  found  using  the  Fourier  transform.  The  signal  as  it 
appears  in  the  frequency-domain  is  displayed  in  Figure  2.  From  the  spectrum,  it  is 
difficult  to  identify  the  number  of  signals,  and  extremely  difficult  to  determine  their 
nature.  Only  an  experienced  signal  analyzer  may  be  able  to  determine  these.  A 
time-frequency  representation  makes  signal  analysis  almost  instantaneous.  The 
Margenau-Hill  spectrogram  (see  Section  3.2)  of  the  signal  is  displayed  in  Figure  3. 
From  the  time-frequency  representation  it  is  immediately  evident  that  there  are  three 
chirp  signals.  The  nature  of  the  chirps  is  easily  determined  as  well.  All  three  chirps 
have  linear  frequency  modulation.  Two  chirps  have  decreasing  frequency,  one  from 
0.45  to  0.05  Hz  and  the  other  from  0.35  to  -0.05  Hz.  The  third  chirp  has  increasing 
frequency  from  -0.5  to  0.5  Flz.  The  time-frequency  representation  can  easily  be 
interpreted  by  someone  with  little  signal  analyzing  experience. 
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Spectra!  Analysis 


Figure  2:  The  spectrum  of  the  multi-component  signal 


Margenau-Hill  Spectrogram  Distribution 
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Figure  3:  The  Margenau-Hill  spectrogram  of  the  multi-component  signal 
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2.2  Extracting  a  Signal  from  Noise 


Due  to  the  non-ideal  conditions  in  collecting  real  data,  signals  will  usually  be  noisy. 
Detecting  the  signal  through  the  noise  is  a  very  important  procedure.  Figure  4  shows  a 
signal  hidden  among  complex  white  Gaussian  noise.  The  signal  was  created  by  adding 
noise  to  a  generated  chirp  signal. 


Real  Part  of  Signal  In  Time 


Figure  4:  Time  series  of  a  signal  in  white  Gaussian  noise 


From  the  time  series  it  is  difficult  to  detect  the  signal  and  nearly  impossible  to 
determine  its  nature.  The  spectrum  of  the  signal,  shown  in  Figure  5,  is  difficult  to  read. 
The  target  signal  is  indistinguishable  from  the  noise. 

Looking  at  a  time-frequency  representation  of  the  signal  in  noise,  it  is  clear  that  the 
signal  is  a  single  chirp.  In  Figure  6,  the  signal  can  be  seen  arching  across  the 
time-frequency  plane.  The  spectrum  of  the  signal  does  not  provide  an  accurate  method 
of  analyzing  the  signal.  In  order  to  properly  identify  the  noisy  target  signal,  a 
time-frequency  representation  is  needed. 

Once  the  signal  is  identified,  the  noise  can  be  filtered  out.  In  this  case,  the  transform  is 
run  through  a  high  pass  filter.  The  filtered  signal  can  then  be  time-averaged  to  convert 
the  signal  to  its  one-dimensional  frequency-domain  representation.  To  convert  this 
frequency-domain  signal  back  to  the  time-domain,  the  inverse  Fourier  transform  is 
taken.  The  filtered  time-frequency  representation  and  the  recovered  signal  are  shown  in 
Figure  7.  The  recovered  signal  is  black  and  the  original  signal  before  noise  was  added 
is  grey.  Using  more  accurately  matched  filters  will  increase  the  quality  of  the  recovered 
signal. 
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Spectral  Analysis 


Figure  5:  The  spectrum  of  a  signal  in  white  Gaussian  noise 


Figure  6:  A  time-frequency  representation  of  a  signal  in  white  Gaussian  noise 
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Filtered  Gabor  Representation 
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Figure  7:  Filtered  time-frequency  transform  and  time  signals 
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One  of  the  applications  of  the  time-frequency  analysis  is  the  detection  of 
noise-corrupted  signals.  In  general,  random  noise  tends  to  spread  evenly  across  the 
time  and  frequency  domains.  However,  the  signal  usually  concentrates  in  a  relatively 
short  period  or  a  narrow  frequency  band.  If  one  converts  the  noise-corrupted  signal  to 
the  joint  time-frequency  domain,  one  can  substantially  improve  the  local  or  regional 
Signal-to-Noise  (SNR). 

2.3  Signal  Analysis 

This  section  will  show  the  usefulness  of  time-frequency  analysis  in  analyzing  real 
aircraft  data.  The  data  had  been  collected  by  Pulse  Doppler  High  Frequency  Surface 
Wave  Radar,  which  used  a  10-element  linear  receiving  antenna  array.  The  data  had 
been  collected  during  the  known  presence  of  a  target.  The  radar  carrier  frequency  is 
5.672  MHz  and  the  pulse  repetition  frequency  is  9.17762  Hz.  The  Fourier  spectrum  of 
this  signal  is  shown  in  Figure  8.  The  central  portion  of  the  spectrum,  which  comprises 
two  strong  spectral  spikes,  is  occupied  by  the  sea-clutter.  From  the  spectrum,  it  appears 
that  the  target  is  present  between  0.3  and  0.4  Hz.  The  spread  of  its  spectrum  tells  us 
very  little  about  the  target’s  behavior. 

The  time-frequency  representation  immediately  provides  more  information.  Figure  9 
shows  the  Gabor  representation  of  the  signal.  From  Figure  9,  the  target  signal  can  be 
seen  increasing  from  slightly  below  0.3  Hz  to  0.4  Hz.  The  increase  is  non-uniform  and 
occurs  between  time  bins  1 50  and  225.  This  signifies  that  the  target  is  accelerating 
moderately. 

The  original  target  signal  can  easily  be  recovered  from  the  time-frequency 
representation.  First  the  transform  is  filtered  to  obtain  only  the  target  signal.  The 
filtered  transform  is  then  averaged  over  time  to  obtain  the  target  signal  in  the 
frequency-domain.  The  frequency  of  the  signal  is  then  converted  to  the  time-domain  by 
an  inverse  Fourier  transform.  The  filtered  transform  can  be  seen  in  Figure  10.  The 
recovered  signal  compared  to  the  original  sampled  data  can  be  seen  in  Figure  1 1 .  The 
recovered  signal  is  black  and  the  original  signal  is  grey. 

Through  the  use  of  time-frequency  representation,  the  analysis  of  the  signal  is  greatly 
simplified.  The  image  of  the  transform  allows  for  quick,  visual  analysis  of  the  target 
signal.  Through  the  use  of  a  filter,  the  original  target  signal  can  be  easily  recovered  for 
further  signature  and  detection  analysis.  When  the  SNR  is  very  low,  as  with  many  radar 
signals,  the  joint  time-frequency  analysis  might  offer  the  only  opportunity  to  detect  the 
signal  of  interest. 
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Gabor  Representation 
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Figure  10:  Filtered  Gabor  representation 


Figure  11:  Recovered  target  signal 
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3.  Comparison  of  Time-Frequency  Distributions 


There  are  many  time-frequency  distributions  that  have  been  developed.  The  transforms 
we  will  be  examining  include  linear,  bilinear,  adaptive,  and  averaged  transforms.  All 
the  time-frequency  representations  can  be  obtained  from 


(1) 


Cx{t,v;f)  =  j  J  J  exp-j27r^4'  f(£,T)x(s+T/2)x*(s-T/2)exp  j2nuT  d^dsdr 


where  x  is  the  signal,  t  is  the  time,  v  is  the  frequeny  and  r  is  the  running  time.  /(£,  r) 
is  a  two-dimensional  function  called  the  parameterization  function  or  the  kernel  defined 
along  the  £-axis  and  r-axis.  This  class  of  distribution  is  known  as  the  Cohen’s  class 
[5-6,1 1-12],  which  can  be  written  as 


(2)  Cx(t,v',U)  =  j  Jn(s-t,t-v)Wx{8,t)d8dt 
where 

(3)  II(t,i/)=  J  J  f(Z,T)  exp-j2v(uT+W  dtdv 

is  a  two-dimensional  Fourier  transform  of  the  kernel  function  /  and  Wx(s,  £)  is  the 
Wigner-Ville  distribution.  This  class  is  significantly  important  since  it  includes  a  large 
number  of  the  existing  time-frequency  energy  distributions.  Moreover,  this  general 
class  can  be  used  to  facilitate  the  design  of  new  time-frequency  transform.  There  are  a 
number  of  alternative  ways  of  writing  the  general  class  of  time-frequency  distributions 
that  are  convenient  and  add  considerable  physical  insight  into  the  process.  The 
properties  of  a  distribution  are  reflected  by  simple  constraints  on  the  kernel,  by 
examining  the  kernel,  one  can  readily  ascertain  the  properties  of  the  distribution.  This 
allows  one  to  pick  and  choose  those  kernels  that  produce  distributions  with  prescribed, 
desirable  properties.  Moreover,  the  Wigner-Ville  distribution  is  the  element  of  the 
Cohen’s  class  for  which  the  function  II  is  a  double  Dirac:  II(t,  v)  =  6{t)8(y),  i.e., 

/(£,  r)  =  1.  For  the  Wigner-Ville  distribution  the  kernel  is  one. 

In  this  section  we  will  begin  with  background  information  on  the  twelve  transforms  that 
are  being  evaluated.  This  background  provides  a  short  description  of  each 
time-frequency  transform.  Detailed  descriptions,  properties  and  derivations  of  the 
transforms  can  be  found  from  the  references,  particularly  from  [5-7,12-14]. 

The  performance  of  the  twelve  transforms  will  then  be  compared  through  five 
experimental  radar  data  sets.  These  sets  will  range  from  returns  with  constant 
frequency  to  returns  with  non-uniform  frequency  modulation.  Each  transform  will  then 
be  evaluated  for  its  effectiveness  in  various  situations. 
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3.1  Gabor  Representation 

The  Gabor  representation  performs  an  expansion  on  a  time  signal.  The  expansion  is 
represented  as  a  finite  set  of  time-shifted  and  frequency  modulated  signals.  These 
signals  are  localized  in  a  time-frequency  plane.  The  Gabor  expansion  can  be  calculated 
with: 


(4)  Gx  [n,  m;  h]  —  x[fc]  h*  [k  —  n]  exp  2j,7rTlfc 

k=—oo 

where  x  is  the  signal,  h  is  the  synthesis  window,  m  is  the  number  of  Gabor  coefficients 
in  frequency  and  n  is  the  number  of  Gabor  coefficients  in  time.  The  Gabor 
representation  is  a  linear  time-frequency  representation  [15-16], 

3.2  Margenau-Hill  Spectrogram 

The  Margenau-Hill  spectrogram  is  defined  as: 


(5)  MHSx(t,  v)  =  U{K^Fx(t,  g)Fx(t,  v;  h)} 

where  Kgh  =  f  h(u)g*(u)du  and  Fx(t,  v\  g)  is  the  short-time  Fourier  transform  of  x. 
h  and  g  are  the  analysis  windows.  The  Margenau-Hill  spectrogram  is  a  bilinear 
time-frequency  representation  [17-18], 

3.3  Born-Jordan  Distribution 

The  Born-Jordan  distribution  is  a  shift-invariant,  kernel  smoothed  Wigner-Ville 
distribution.  The  kernel  function  defined  for  equation  (3)  is 


(6) 


m,r) 


sin(7r£r) 
7 r£r 


which  defines  the  Born-Jordan  distribution  as: 


/•oo  1  /-t+M/2 

(7)  BJx(t,u)=  /  — t  /  x(s  +  t/2)x*(s  —  r/2)ds.  exp~J  7ri/T  dr 

J— oo  r  I  Jt— |r|/2 

where  t  is  the  time,  u  is  the  frequency,  x  is  the  signal,  and  r  is  the  running  time.  The 
conjugation  of  x  is  x*.  The  Born-Jordan  distribution  is  a  bilinear  time-frequency 
representation  [11], 
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3.4  Binomial  Distribution 


The  binomial  distribution  is  a  shift-invariant,  reduced  interference  distribution  with  a 
kernel  based  on  the  binomial  coefficients.  The  binomial  distribution  is  defined  as: 


°°  ,  M 

RIBNx(t,  v)  =  Y,  73  E 

T=  —  OO  V=-\ T 

(8) 


22|t|+! 


2M  +  1 

r  +  v  +  1 


x[t+v+T]x*[t+v-T]  exp_j47r"T 


The  binomial  distribution  is  a  bilinear  time-frequency  representation  [19], 


3.5  Choi-Williams  Distribution 

The  Choi-Williams  distribution  was  developed  as  a  method  of  finding  the  Wigner-Ville 
distribution  with  the  minimum  amount  of  cross-term  interference.  The  Choi-Williams 
distribution  is  a  shift-invariant  transform.  Essentially,  the  distribution  is  a  smoothed 
version  of  the  Wigner-Ville  distribution  through  a  kernel  function  defined  by 


(9) 


/(£>t)  =  exp 


2<t2 


The  Choi-Williams  distribution  is  then  defined  as 


CWx(t,v)  =  yjl  J  J°°  y^w-2^s~t)2/T2  x{s+^)x*(s-^)exp-j2™r  dsdr 
(10) 

The  smoothing  of  the  distribution  is  controlled  by  the  constant  a.  As  a  — >  oo,  the 
Choi-Williams  distribution  will  simplify  converge  to  the  Wigner-Ville  distribution,  as 
the  kernel  goes  to  1 .  The  cross-terms  are  reduced  with  smaller  values  of  o.  The 
Choi-Williams  distribution  is  a  bilinear  time- frequency  representation  [20]. 


3.6  Smoothed  Wigner-Ville  Distribution 

The  smoothed  pseudo  Wigner-Ville  distribution  is  given  as 


SPWVx(t,u;g,  h) 

(11) 


g(s-t)x(s-\-T /2)x* (s-t j2)ds exp  ■?27n'T dr 
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where  h  is  the  frequency  smoothing  window  and  g  is  the  time  smoothing  window.  The 
smoothed  Wigner-Ville  distribution  considered  in  this  analysis  is  Gaussian  “smoothed". 
Smoothing  the  Wigner-Ville  distribution  causes  a  loss  of  phase  information  but  it  can 
eliminate  to  a  certain  degree  the  cross-term  interference.  The  smoothed  Wigner-Ville 
distribution  is  a  bilinear  time-frequency  representation  [13-14], 

Spectrogram  Distribution 

The  spectrogram  distribution  is  the  square  of  the  short-time  Fourier  transform.  It  is  also 
known  as  the  energy  density  spectrum.  An  advantage  of  the  spectrogram  distribution  is 
that  it  yields  all  positive  values.  The  spectrogram  is  given  as 


(12) 


SPx(t,  v)  = 


/*(«)*>- ()e* 


where  h  is  the  frequency  smoothing  window  and  h  being  normalized  so  as  to  be  of  unit 
energy.  The  spectrogram  distribution  is  a  bilinear  time-frequency  representation 
[6,13-14]. 


Zhao-Atlas-Marks  Distribution 

If  we  smooth  the  Bom-Jordan  distribution  (section  3.3)  along  the  frequency  axis,  we 
obtain  the  Zaao-Atlas-Marks  distribution.  The  distribution  is  defined  as 


(13) 


ZAMx(t,v) 


ft+  M/2 

h(r )  /  x(s  +  t/2)x*(s  —  r/2)ds 

Jt-\r\/2 


exp  ■?27r"r  dr 


It  is  also  known  as  the  Cone-Shaped  kernel  distribution.  The  Zhao-Atlas-Marks 
distribution  is  a  bilinear  time-frequency  distribution  [21], 


Reassigned  Gabor  Representation 

The  reassigned  time-frequency  representation  performs  an  averaging  method  upon  their 
base  transform  [13],  The  average  values  are  shifted  to  the  center  of  gravity  of  the 
energy  contributions.  Through  the  averaging  and  shifting,  bilinearity  of  the  transform  is 
lost.  Time  and  frequency  shift  invariance  is  preserved  and  the  energy  reassignment 
obeys  the  conservation  of  energy  requirement.  The  reassigned  Gabor  representation  at 
point  (t' ,  v)  is  given  by  the  following  expression 


(14)  Sx\t  ,v  \h)  =  J  J  Sx(t,u;  h)6(t  —  i(x;t,v))8(u  —  v{x\  t ,  v))dtdv 
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where  Sx  is  the  spectrogram  time-frequency  distribution.  The  analysis  window  h  used 
in  this  spectrogram  is  a  Gassian  window.  The  reassignment  method  moves  each  value 
of  the  spectrogram  computed  at  any  point  ( t ,  v)  to  another  point  (t' .  v),  which  is  the 
center  of  gravity  of  the  signal  energy  distribution  around  ( t ,  v).  t(x;  t,  v)  and  v(x\  t,  v) 
are  defined  as; 


(15) 


and 


(16) 


t{x\  t,v)  =  t  —  3? 


[  Fx(t,u-,Th)F*(t,u;  ft)) 
\  \Fx(t,  w,  h)|2  J 


v(x;  t,  v)  =  V  +  Q 


[Fx{t,v-Dh)F*{t,v-h)\ 

\  2'K\Fx{t,v-h)\2  J 


As  well,  Th(t )  =  t  hit)  and  Dh(t)  —  %{t).  The  windows  Th  and  Dh  defined  above  are 
colinear  in  this  case. 


3.10  Reassigned  Morlet  Scalogram 

The  Morlet  scalogram  is  the  square  of  the  wavelet  transform  using  Morlet  wavelets 
[13],  The  wavelet  transform  uses  a  windowing  process  different  to  the  sliding  window 
used  in  the  short-time  Fourier  transform.  The  wavelet  transform  uses  narrow  windows 
at  high  frequencies  and  wide  windows  at  low  frequencies.  The  reassigned  Morlet 
scalogram  can  be  calculated  as: 


(17)  Sdr\t,a;h)  -  J  j  a'2SCx{t,a;h)6(t 


we/'  w  dtda 

t{x\t,a))b{a  -a(x;t,a)) — — 


where 


(18) 


and 


(19) 


t(x 


.  t  \  f  {Tx(t,a;Th)T*(t,a; 


h) 


i>(x-,t,  a) 


vo  _  ,  o  /  TxiPi  a'i  F>h)Tx  ( t ,  a;  h)  ) 

a(x-,t,a)  a  ^  \  2'rra\Tx(t,a-,h)\‘2  j 


with  Th(t )  =  t  h(t )  and  Dh(t)  =  %{t),  where  h(t )  is  a  Gassian  window.  SCx{t,  a;  h ) 
is  the  scalogram  and  Tx(t,  a;  h )  is  the  wavelet  transform: 


DRDC  Ottawa  TM  2003-090 


15 


(20) 


SCx(t,a;h )  =  \Tx(t,  a;  h)\2  = 


J  x(s)h*(^-)di 


3.11  Reassigned  Spectrogram  Distribution 

The  reassigned  spectrogram  distribution  performs  an  averaging  and  shifting  method 
upon  the  spectrogram  distribution  [13],  The  reassigned  spectrogram  is  defined  as: 


(21)  Sir\t,v';h)=  J  J  Sx(t,v,ti)6(t  -  t(x;t,t/))6(u' 


i>{x\  t,  v))dtdv 


where 


(22) 


and 


t(x;t,v)  =  t  -  5? 


\Fx{t,v-Th)F*(t,u-h)\ 
X  I  Fx(t,ir,h)\*  J 


(23) 


i>(x;  t,u)  =  v  +  9 


f  Fx{t,u\Dh)F*{t,v,h)\ 
X  2n\Fx(t,  u;  h)\2  / 


Sx(t,  v\  h)  is  the  spectrogram  distribution  and  Th(t)  =  t  h(t )  and  D/t(t)  =  fw- 


3.12  Adaptive  Energy  Distribution 

Adaptive  energy  distribution  (AED)  is  a  window-matched  spectrogram,  in  which  the 
window  matching  is  done  with  the  use  of  generalized  instantaneous  parameters  (GIP’s). 
The  GIP’s  that  are  used  for  Gaussian  window  matching  are  the  instantaneous  time 
variance  (ITV),  the  instantaneous  frequency  variance  (EFV),  and  the  instantaneous 
time-frequency  (ITF)  covariance  [22], 


(24) 


ITVgh(t,is)  =  ~3t 


d2Sx(t,v)  1 
dv2  Sx{t,v) 


(dSx{t,v)  1  \2 

V  df  Sx(t,  v) ) 


(25) 


IFVzh(t,u)  =  -~R 


d2Sx(t,  u)  1 
dt2  Sx(t,v) 


(dsx{t,v)  1  \2 

V  dt  Sx{t,v)J 
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TTF  ( t  _ _ 1  ^\d2sx{t,v)  i  dsx{t,v)dsx{t,v)  f  1  \2 

zU  ’  ;  8tt2  Sx(t,i/)  01/  dt  \Sx(t,  v) ) 


where  Sx(t,  v )  specifies  the  spectrogram  of  the  signal  x(t)  formed  using  the  window 
h(t).  The  subscript  zh  in  equations  (24)-(26)  denotes  that  the  signal  x(t)  has  been 
multiplied  by  the  window  h(t).  The  parameters  for  an  optimal  window  are  found  by 
matching  the  GIP' s  of  the  signal-window  spectrogram  to  those  of  the  window-window 
spectrogram.  The  optimal  window  parameters  are  updated  iteratively  until  the  desired 
level  of  matching.  At  the  nth  iteration  they  are  calculated  from 


an+i(t,u) 


(8n)3(ITFzhn(t,u)flFVzhn(t,v) 
(8n)2(ITFzhn(t,  is))2  +  1 


Pn+l(t>v)  = 


(8tt  )2ITFzhn(t,v)IFVzhn(t,v) 
(8n)2(ITFzhn(t,  u))2  +  1 


or  if  ITVzhn(t,  v )  >  IFVzhn(t,  v)  then 


an+i(t,v)  = 


8t t  ■  ITVzhn(t,u) 


Pn+l(Pi  v)  = 


(8tt  )2ITVzhn(t,u)ITFzhn(t,u) 


where  a  and  /3  are  the  parameters  of  the  Gaussian  window  h(t)  =  exp  n(a 
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4.  Comparison  of  Transforms _ 

4.1  Comparison  Signals 

The  signals  in  the  following  comparisons  are  experimental  aircraft  data.  These  signals 
were  collected  by  a  High  Frequency  Pulse  Doppler  Surface  Wave  Radar,  which  used  a 
10-element  linear  receiving  antenna  array.  The  data  was  collected  with  a  target  present. 
The  radar  carrier  frequency  is  5.672  MHz  and  the  pulse  repetition  frequency  is  9.17762 
Hz.  There  are  5  trials  and  each  trial  corresponds  to  a  block  of  256  pulses.  The  detailed 
description  of  the  radar  is  given  by  [23,24]. 


Signal  1  Signal  2 


Figure  12:  Spectrums  of  the  five  comparison  signals 


The  Fourier  spectrums  of  the  five  signals  can  be  seen  in  Figure  12.  In  general,  the 
background  against  which  targets  must  be  detected  by  a  radar  is  clutter  and  noise 
(thermal  and  atmospheric).  In  the  case  of  High  Frequency  Surface  Wave  Radar 
(HFSWR)  used  over  the  ocean,  the  clutter  is  due  to  Bragg  scattering  from  the  surface  of 
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the  ocean  and  is  known  as  sea-clutter  [23,24],  The  sea-clutter  components  are 
concentrated  around  zero-Doppler  with  two  large  spectral  spikes  around  zero-Doppler 
known  as  the  Bragg  lines.  The  spectral  peak  on  the  left  hand  side  of  Figure  12  is  called 
the  receding  Bragg  line.  The  one  on  the  right  hand  side  is  called  the  advancing  Bragg 
line.  The  Bragg  lines  appear  as  two  veiy  large  targets.  Nomally  high  speed  aircraft 
appears  in  the  flat  portion  of  the  Doppler  spectrum.  However,  when  the  aircraft  travels 
in  a  direction  tangential  to  the  radar  beam,  as  in  the  cases  of  signals  2  and  3  in  Figure 
12,  the  aircraft  target  can  appear  in  the  velocity  spectrum  close  to  the  zero  value. 


In  Signal  1,  the  target  is  seen  in  the  spectrum.  The  target  is  a  spike  at  around  -0.25  Hz. 
In  Signal  2,  the  target  is  difficult  to  detect.  It  is  near  the  clutter,  but  it  is  spread  and 
difficult  to  distinguish  from  the  clutter.  In  Signal  3,  the  target  is  spread  widely  from 
approximately  -0.3  to  0.2  Hz.  Whether  the  frequency  has  been  increasing,  decreasing, 
or  oscillating  is  undeterminable.  In  Signal  4,  the  target  is  spread  between  -0.4  to  -0.25 
Hz.  Whether  the  frequency  has  been  increasing,  decreasing,  or  oscillating  is  again 
undeterminable.  In  Signal  5,  the  target  is  spread  from  0.35  Hz  to  the  clutter.  Whether 
the  frequency  has  been  increasing,  decreasing,  or  oscillating  is  undeterminable. 


These  five  signals  were  used  to  test  the  performance  of  12  time-frequency  transforms. 
The  parameters,  such  as  smoothing  length  or  kernel  spread,  of  each  time-frequency 
representation  were  chosen  to  provide  the  best  resolution  of  the  target. 


4.2  Constant  Frequency 


The  first  signal  is  a  chirp  with  constant  frequency  modulation.  This  signal  tests  the 
ability  of  the  transforms  to  resolve  constant  frequency  modulation.  The  transforms  are 
shown  in  Figure  13. 


The  target  signal  is  perfectly  resolved  with  the  reassigned  Gabor  representation  and  the 
reassigned  Morlet  scalogram.  The  target  signal  is  nearly  perfectly  resolved  in  time  and 
frequency  with  all  of  the  other  transforms.  Of  these,  the  Choi-Williams  distribution  and 
the  binomial  distribution  provide  the  sharpest  response.  The  reassigned  spectrogram’s 
image  of  the  target  is  very  grainy.  The  adaptive  energy  distribution  and  the  Zhao-Atlas- 
Marks  distribution  are  slightly  grainy  when  compared  with  the  other  transforms. 


The  smoothed  Wigner-Ville  distribution  is  the  “noisiest"  transform.  The  signal  is  very 
faint  compared  to  the  clutter  and  the  interference  cross-term.  The  Choi-Williams 
distribution  is  the  second  noisiest  transform.  While  the  interference  cross-term  is 
reduced,  it  is  still  quite  visible.  The  binomial  distribution  and  the  Bom-Jordan 
distribution  have  small  smears  between  the  signal  and  the  clutter  as  the  only  remnants 
of  the  cross-term.  The  remaining  transforms  are  relatively  noiseless. 
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4.3  Non-Uniform  Chirp  near  Clutter 


The  second  signal  is  a  chirp  with  non-uniform  frequency  modulation  at  normalized 
frequencies  close  to  0  Hz.  This  signal  tests  the  ability  of  the  transforms  to  detect  a 
target  signal  close  to  the  clutter.  The  transforms  are  shown  in  Figure  14. 

The  Born- Jordan  distribution  and  the  Zhao- Atlas-Marks  distribution  are  unable  to 
clearly  display  the  signal  close  to  the  clutter.  The  target  signal  is  broken  up  and  joined 
with  the  clutter  in  various  places. 

The  binomial  distribution,  the  Choi- Williams  distribution,  and  the  smoothed  Wigner- 
Ville  distribution  display  the  target  signal,  but  it  is  nearly  indistinguishable  from  the 
clutter. 

The  reassigned  transforms  provide  perfect  time  and  frequency  resolution.  The  Gabor 
representation,  the  spectrogram,  and  the  adaptive  energy  distribution  provide  nearly 
perfectly  resolved  signals  in  time  and  frequency,  which  are  easily  distinguishable  from 
the  clutter.  The  Margenau-Hill  spectrogram  provides  very  good  frequency  resolution, 
but  its  time  resolution  is  not  as  good  as  the  other  transforms.  The  reassigned 
spectrogram  and  the  adaptive  energy  distribution  are  grainy  compared  to  the  other 
transforms. 


4.4  Chirp  Through  Clutter 

The  third  signal  contains  a  chirp  with  approximately  linear  frequency  modulation.  This 
chirp  passes  directly  “through"  the  clutter.  This  signal  tests  the  ability  of  the  transform 
to  resolve  weak  signals  that  pass  through  strong  noisy  clutter.  The  transforms  are 
shown  in  Figure  15. 

The  target  signal  is  visible  within  all  of  the  transforms.  With  the  Zhao-Atlas  Marks 
distribution,  the  target  is  seen  only  as  a  disruption  in  the  clutter.  For  the  Bom-Jordan 
distribution,  the  target  signal  is  broken  up  and  smeared  with  the  clutter. 

The  reassigned  Gabor  representation  and  the  reassigned  Morlet  scalogram  perfectly 
resolve  the  target  signal  in  time  and  frequency.  The  clutter  has  minimal  effect  on 
displaying  the  target  signal.  The  smoothed  Wigner-Ville,  the  binomial  distribution,  the 
Choi-Williams  distribution,  and  the  reassigned  spectrogram  also  provide  excellent 
signal  resolution  in  time  and  frequency.  The  cross-term  interference  is  significant  with 
the  binomial  distribution,  Choi-Williams  distribution  and  the  smoothed  Wigner-Ville 
distribution.  Some  parts  of  the  target  signal  have  been  slightly  masked  by  the 
interference.  The  remaining  transforms  provide  decent  frequency  resolution.  The  target 
is  smeared  in  the  frequency-domain  for  these  transforms,  but  it  is  visible  and  does  not 
suffer  from  interference.  The  Gabor  representation  and  the  spectrogram  distribution 
provide  the  best  image  of  the  target  signal  of  these  remaining  transforms. 
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4.5  Changing  Signal 


The  fourth  signal  is  a  chirp  that  changes  from  a  linear  frequency  modulation  to  a 
constant  modulation.  This  signal  tests  the  ability  to  display  a  target  signal  that  changes 
its  characteristic  behavior.  The  transforms  are  shown  in  Figure  16. 

The  target  return  is  perfectly  resolved  in  time  and  frequency  with  the  reassigned  Gabor 
representation  and  the  reassigned  Morlet  scalogram.  The  constant  portion  of  the  signal 
is  faint,  but  is  still  visible  for  both  transforms.  The  smoothed  Wigner-Ville  distribution, 
the  binomial  distribution,  the  Choi- Williams  distribution  and  the  reassigned 
spectrogram  distribution  provide  a  nearly  perfectly  resolved  target  signal  in  time  and 
frequency.  The  constant  portion  of  the  signal  is  faint  with  the  reassigned  spectrogram. 
The  Choi-  Williams  distribution  and  the  smoothed  Wigner-Ville  distribution  have 
significant  cross-term  interference.  The  binomial  distribution  has  a  small  cross-term 
remaining. 

The  Zhao-Atlas-Marks  distribution  does  not  resolve  the  target  clearly.  The  target  is 
smeared  and  grainy.  The  clutter  with  this  distribution  appears  much  larger  than  with  the 
other  distributions. 

The  remaining  transforms  are  a  mix  of  decent  and  good  resolution.  For  the  linear 
frequency  modulation,  the  target  signal  is  well  resolved  in  time  and  frequency.  The 
constant  portion  is  much  sharper.  The  Gabor  representation,  the  Margenau-Hill 
spectrogram  and  the  adaptive  energy  distribution  provide  the  best  images  of  the  target 
of  these  remaining  transforms.  The  Bom-Jordan  distribution  has  cross-term 
interference  approximately  equal  in  magnitude  to  the  target  signal. 


4.6  Small  Frequency  Modulation 

The  final  signal  contains  a  target  signal  with  nearly  linear  frequency  modulation.  The 
modulation  begins  by  increasing  non-uniformly.  It  then  remains  constant  for  a  period, 
and  then  increases  again.  This  signal  tests  the  ability  to  detect  minor  changes  in  the 
behavior  of  a  signal.  The  transforms  are  shown  in  Figure  17. 

The  reassigned  Gabor  representation  and  the  reassigned  Morlet  scalogram  display  a 
perfectly  resolved  target  signal  in  both  time  and  frequency.  All  three  portions  of  the 
target  signal  are  clearly  visible.  The  smoothed  Wigner-Ville,  the  binomial,  and  the 
Choi-Williams  distributions  provide  nearly  perfectly  resolved  target  signals  in  time  and 
frequency.  A  strong  interference  cross-term  is  present  with  the  smoothed  Wigner-Ville 
distribution.  A  smeared  cross-term  is  still  visible  with  the  binomial  distribution.  A 
strong,  smeared  cross-term  is  visible  with  the  Choi-Williams  distribution.  These 
distributions  show  the  three  portions  of  the  signal  clearly. 

The  Zhao-Atlas-Marks  distribution  displays  a  grainy  target  with  poor  frequency 
resolution.  The  clutter  appears  strongest  with  this  transform. 
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The  remaining  transforms  provide  similar  resolution  of  the  target  signal.  While  the 
three  portions  of  the  signal  are  visible  in  all  of  the  transforms,  certain  transforms  do  not 
display  them  very  well.  The  curves  of  the  target  signal  are  less  pronounced  with  the 
Gabor  representation,  the  Margenau-Hill  spectrogram,  and  the  adaptive  energy 
distribution,  which  make  the  target  return  appear  to  be  a  linear  chirp.  The  Bom- Jordan 
distribution  displays  the  nature  of  the  signal  well,  but  it  is  grainy  and  smeared  along  the 
frequency  axis.  The  spectrogram  and  the  reassigned  spectrogram  provide  good  time 
and  frequency  resolution,  but  the  reassigned  spectrogram  is  very  grainy. 
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Figure  13:  Time-frequency  representations  of  a  signal  with  constant  frequency. 
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Figure  14:  Time-frequency  representations  of  a  signal  near  clutter. 
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Figure  15:  Time-frequency  representations  of  a  signal  through  clutter. 
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Figure  16:  Time-frequency  representations  of  the  changing  signal. 
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Figure  17:  Time-frequency  representation  of  a  signal  with  small  frequency  modulation. 


DRDC  Ottawa  TM  2003-090 


27 


5.  Evaluation  of  Transforms 

5.1  Gabor  Representation 


The  Gabor  representation  provides  excellent  resolution  in  time  and  frequency  for  all 
signals  except  those  with  a  large  chirp.  Signals  close  to  the  clutter  are  not  affected  by 
the  clutter.  No  cross-term  interference  is  present  in  the  Gabor  representation.  The 
Gabor  representation  has  difficulty  in  sharply  resolving  chirp  signals.  The  signal  is 
smeared  for  large  frequency  changes.  For  cases  where  frequency  changes  are  small,  the 
Gabor  representation  is  a  good  transform. 


5.2  Margenau-Hill  Spectrogram 

The  Margenau-Hill  spectrogram  provides  good  resolution  in  time  and  frequency  for  all 
signals  except  the  chirp  through  the  clutter.  For  signals  close  to  the  clutter,  the 
transform  loose  some  time  resolution.  This  transform  has  no  cross-term  interference. 
For  cases  where  frequency  changes  are  small,  the  Margenau-Hill  spectrogram  is  a  good 
transform. 

5.3  Born-Jordan  Distribution 

The  Born-Jordan  distribution  performs  well  for  signals  with  constant  frequency.  The 
signal  is  warped  and  broken  up  for  signals  close  to  the  clutter.  The  distribution  displays 
smeared  cross-term  interference,  although  it  is  generally  weaker  than  the  signal.  For 
signals  with  frequency  modulation,  the  Born-Jordan  distribution  does  not  provide  a 
smooth  target  image.  This  transform  is  not  ideal  due  to  its  inability  to  resolve  targets 
near  clutter  and  its  poor  resolution  with  chirp  signals. 

5.4  Binomial  Distribution 

The  binomial  distribution  provides  good  resolution  in  time  and  frequency  for  signals  of 
all  types.  The  signal  is  still  visible  close  to  clutter,  although  cross-tenn  interference 
impedes  easy  signal  detection.  Smeared  cross-term  interference  is  present  in  the 
binomial  distribution.  Its  magnitude  ranges  from  weaker  than  the  signal  to  stronger 
than  the  signal  depending  on  the  proximity  of  the  signal  to  the  clutter.  For  cases  where 
the  cross-term  interference  is  not  important  and  the  target  is  not  close  to  the  clutter,  the 
binomial  distribution  is  a  good  transform. 

5.5  Choi-Williams  Distribution 

The  Choi-Williams  distribution  provides  excellent  resolution  in  time  and  frequency  for 
signals  of  all  types.  Signals  close  to  clutter  are  difficult  to  detect  due  to  interference. 
Smeared  cross-term  interference  is  present  in  the  Choi-Williams  distribution,  with  a 
magnitude  ranging  from  slightly  less  to  greater  than  the  magnitude  of  the  signal.  For 
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cases  where  the  cross-term  interference  is  not  important  and  the  target  is  not  close  to 
the  clutter,  the  Choi- Williams  distribution  is  a  good  transform. 


5.6  Smoothed  Wigner-Ville  Distribution 

The  smoothed  Wigner-Ville  distribution  provides  excellent  time  and  frequency 
resolution  in  all  cases.  Signals  close  to  the  clutter  are  visible  but  cross-term  interference 
makes  detection  difficult.  The  cross-term  interference  is  prominent  in  all  signals,  with 
large  cross-terms  with  magnitudes  greater  than  the  signal.  For  signals  that  are  not  close 
to  the  clutter  and  cross-terms  are  not  important,  the  smoothed  Wigner-Ville  distribution 
is  an  excellent  transform. 

5.7  Spectrogram  Distribution 

The  spectrogram  distribution  provides  good  resolution  in  cases  of  constant  frequency 
and  small  frequency  changes.  Signals  close  to  the  clutter  show  no  interference.  There  is 
no  cross-term  interference  present  in  this  transform.  For  signals  that  have  small 
frequency  changes,  the  spectrogram  distribution  is  a  good  transform. 

5.8  Zhao-Atlas-Marks  Distribution 

The  Zhao-Atlas-Marks  distribution  provides  good  resolution  in  the  case  of  a  constant 
frequency.  Signals  close  to  the  clutter  are  broken  up  and  combined  with  the  clutter  in 
some  places.  Smeared  cross-term  interference  is  visible  with  this  transform.  For  signals 
with  changing  frequency,  the  Zhao-Atlas-Marks  distribution  displays  ghosting  around 
the  signal.  The  ghosting  increases  when  the  rate  of  frequency  change  increases.  The 
Zhao-Atlas-Marks  is  not  an  ideal  transform  due  to  its  poor  performance  with  signals 
close  to  the  clutter  and  signals  with  changing  frequency. 

5.9  Reassigned  Gabor  Representation 

The  reassigned  Gabor  representation  provides  excellent  time  and  frequency  resolution 
in  all  cases.  The  signal  crossing  the  clutter  is  not  affected  by  the  clutter.  Signals  close  to 
the  clutter  are  broken  up  and  joined  with  the  clutter  in  some  places.  There  is  no 
cross-term  interference  in  this  transform.  For  signals  that  are  not  close  to  the  clutter,  the 
reassigned  Gabor  representation  is  an  excellent  transform. 

5.10  Reassigned  Morlet  Scalogram 

The  reassigned  Morlet  scalogram  provides  excellent  time  and  frequency  resolution  in 
all  cases.  Signals  close  to  the  clutter  are  not  affected  by  the  clutter.  There  is  no 
cross-term  interference  in  this  transform.  The  reassigned  Morlet  scalogram  is  an 
excellent  transform. 
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5.11  Reassigned  Spectrogram  Distribution 

The  reassigned  spectrogram  distribution  provides  very  good  time  and  frequency 
resolution  in  all  cases.  Signals  close  to  the  clutter  are  not  affected  by  the  clutter.  There 
is  no  cross-term  interference  in  this  transform.  The  image  produced  by  the  reassigned 
spectrogram  distribution  is  not  smooth.  The  reassigned  spectrogram  distribution  is  a 
very  good  transform. 

5.12  Adaptive  Energy  Distribution 

The  adaptive  energy  distribution  performs  quite  well  for  signals  with  constant 
frequency  and  small  frequency  changes.  It  avoids  interference  when  the  target  signal  is 
close  to  or  passing  through  clutter.  The  adaptive  energy  distribution  also  does  not 
display  any  cross-term  interference.  The  transform  has  poor  frequency  resolution  for 
signals  with  larger  frequency  changes.  For  cases  with  small  or  no  frequency 
modulation,  the  adaptive  energy  distribution  is  a  good  transform. 
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6.  Overall  Evaluation 


For  the  best  possible  resolution,  the  reassigned  transforms  should  be  used.  The 
averaging  algorithm  guarantees  signals  perfectly  localized  in  frequency.  The  only 
downside  to  these  transforms  is  the  time  involved  for  their  calculation.  Before 
reassignment,  the  base  transform  must  be  performed.  The  reassignment  can  increase 
the  total  time  for  calculation  by  a  factor  of  two  or  more. 

In  cases  where  resolution  is  not  the  main  concern,  the  Gabor  representation,  the 
Margenau-Hill  spectrogram,  and  the  spectrogram  distribution  are  good  transforms  to 
use.  They  are  relatively  quick  compared  to  the  other  transforms.  While  they  do  not 
provide  great  resolution  in  cases  of  large  frequency  modulation,  they  accurately  display 
the  signal. 

The  binomial  distribution,  the  Choi-Williams  distribution,  and  the  smoothed  Wigner- 
Ville  distribution  provide  very  good  resolution  and  good  speed  at  the  expense  of  some 
interference  cross-terms.  For  target  signals  away  from  the  clutter,  the  target  and  the 
cross-term  are  easy  to  identify.  If  the  target  is  close  to  the  clutter,  the  cross-terms  make 
detection  of  the  signal  difficult. 
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7.  Conclusion 


Joint  time-frequency  representations  provide  excellent  ways  of  analyzing 
non-stationary  signals.  Their  ability  to  display  frequencies  that  change  in  time  is  an 
improvement  over  Fourier  analysis  of  signals.  This  study  also  shows  that 
time-frequency  transforms  provide  additional  insight  into  the  analysis,  interpretation, 
and  processing  of  radar  signals  that  are  sometimes  superior  to  what  is  achievable  in  the 
traditional  time  or  frequency  domain  alone. 

With  the  Fourier  spectrum,  multi-component  signals  were  difficult  to  detect,  much  less 
analyze.  With  time-frequency  representations,  the  components  of  a  signal  were 
immediately  visible  due  to  the  ability  to  show  time-varying  frequency. 

Noise  in  signals  greatly  hampered  the  Fourier  spectrum’s  ability  to  detect  and  analyze 
signals.  While  noise  affected  target  detection  with  the  time-frequency  analysis,  its 
effect  was  much  less.  The  noise  could  be  filtered  out  relatively  easily,  which  allowed 
for  the  recovery  of  the  target  signal. 

The  ability  of  time-frequency  analysis  to  analyze  signals  is  demonstrated  using  a  real 
signal  from  HF  radar.  The  process  from  identification  to  recovery  of  the  target  signal 
was  shown. 

Several  time-frequency  transforms  have  been  compared  and  evaluated  using 
experimental  HF  radar  data  for  detecting  maneuvering  air  targets  in  sea-clutter 
environment.  Each  set  of  data  tested  the  ability  of  the  transforms  to  resolve  different 
types  of  signals.  While  each  transform  had  their  own  strengths  and  weaknesses,  several 
transforms  stood  out  for  their  flexibility  and  target  image  resolution. 

From  the  data  used  in  this  report,  the  reassigned  transforms  provided  the  best  signal 
resolution  at  the  expense  of  computational  time.  The  Gabor  representation,  the 
Margenau-Hill  spectrogram  and  spectrogram  distribution  provided  a  quick  alternative 
for  simple  analysis  at  the  expense  of  resolution.  The  binomial  distribution,  the 
Choi-Williams  distribution,  and  the  smoothed  Wigner-Ville  provided  a  good  balance  of 
speed  and  resolution  at  the  expense  of  some  cross-term  interference.  These  results 
demonstrated  that  joint  time-frequency  analysis  techniques  provide  an  effective  method 
of  detecting  maneuvering  air  targets  in  heavily  cluttered  regions. 
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